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ABSTRACT 

This paper describes a practical model for the polarization of Be stars which can 
be used to estimate roughly the physical parameters for optically thin circumstellar 
envelopes from broadband UBVRI photopolarimetry data. Analysis of long-term vari- 
ability in terms of these parameters is a promising approach toward understanding the 
Be phenomenon. An interesting result from fitting the model to observations of eight 
Be stars is that all of them may have geometrically thin disks, with opening half-angles 
on the order of ten degrees or less. This contributes to the growing evidence that most 
Be disks are geometrically thin. 

Subject headings: polarization — scattering — techniques:polarimetric — stars: emission- 
line, Be — stars:winds, outflows — circumstellar matter 

1. Motivation 

The idea that the Ha emission lines of Be stars originate from an equatorially flattened cir- 
cumstellar envelope or disk has now been directly verified by interferometric imaging (Quirrenbach 
et al. 1997). The related explanation of the linear polarization of Be stars as due to scattering of 
the starlight by free electrons in the disk was also confirmed, because the observed position angle 
of the polarization was found to be perpendicular to the direction of elongation of the disk as pro- 
jected onto the plane of the sky. It is this apparent asymmetry of the scattering region that results 
in a net polarization, by defeating the cancellation which would take place if there were identical 
scattering subregions in adjacent quadrants of the projected disk (see Fig. 1). 

A thorough historical review of both observation and theory of polarization in Be stars was 
given by Coyne & McLean (1982). Since that time the major observational developments in the 
field have been increased spectral resolution by spectropolarimetry and extension of the wave- 
length coverage into the ultraviolet by space-based observations (Bjorkman 2000). Beginning with 



Astronomical Institute, University of Amsterdam, Netherlands 



- 2 - 




Fig. 1. — Scattering geometry in the disk of a Be star, illustrating by extremes how the net polarization depends 
on the orientation of the disk on the sky. In the upper sketch the disk is viewed pole-on, with angle of inclination to 
the line of sight i = 0°, while in the lower sketch the disk is viewed equator-on, with i = 90°. Double arrows show the 
direction of vibration of the electric field of starlight scattered in the disk. For the pole-on case cancellation occurs 
between adjacent quadrants of the disk, and the net polarization is zero because the electric fields of the scattered 
light are perpendicular. When the disk is viewed equator-on there is no cancellation, so the observed polarization is 
a maximum. 
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the pioneering work of Poeckert & Marlborough (1978a,b), who calculated the single scattering 
polarization in an NLTE hydrogen envelope, theoretical analysis and modeling of the continuum 
polarization of Be star disks has progressed through increasingly complicated treatments. Fox 
(1991) developed analytic equations to calculate the single scattering polarization for axisymmetric 
envelopes. Bjorkman h Bjorkman (1994) extended the single scattering calculations to include at- 
tenuation and emission within the envelope. Hillier (1994) developed polarization source functions 
using traditional Feutrier methods to integrate the transfer equation for optically thick cases, while 
Wood, Bjorkman, & Bjorkman (1997, hereafter WBB) developed Monte Carlo methods. 

One may wonder, then, if there is anything new to be learned about Be stars from continued 
broadband polarimetry. The answer is, of course, that even though we now have firmer knowledge 
of the disk characteristics, including temperature, density, and geometry, broadband polarization 
monitoring is still a primary source of information on the variability of these quantities over time 
scales of years to decades, which are typical of the most fundamental aspect of the Be phenomenon: 
the unpredictable transition from the normal B phase to the Be phase and back again, associated 
with the formation and later dissipation of the circumstellar envelope. 

This paper attempts to establish a simple and approximate, but nevertheless serviceable ap- 
proach to the modeling and analysis of UBVRI polarimetry data to derive the basic envelope 
parameters. Tracing the observed variability in terms of these parameters still promises to give 
valuable hints about the nature of the unknown processes involved in the Be phenomenon. 



Figure 2 shows the geometry of the circumstellar envelope model adopted for this study. Intro- 
duced by Kruszewski, Gehrels, & Serkowski (1968) to investigate the polarization of red variables, 
it was later used by Brown & McLean (1977) to illustrate their theoretical formulation of the po- 
larization by electron scattering in a Be disk. Its shape may be described as an axially symmetrical 
sector of a sphere, with a wedge-shaped cross section opening outward at half-angle a. 

For small opening angles the spherical sector is a good representation of an equatorial disk, 
with the advantage of being mathematically suited for spherical coordinates. A basic starting point 
is to assume the disk is pure hydrogen with uniform electron temperature T e , extending to infinite 
distance (R e — > oo in Fig. 2) from a central star of radius i?*, with electron number density N e 
given by a radial power law with exponent rj : 



2. A Spherical Sector Envelope Model 




(1) 



Waters, Cote, &; Lamers (1987) derived values of 2.0 < rj < 3.5 for some of the stars also 
studied in this paper using the slope of the infrared continuum from IRAS observations, but the 
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Fig. 2. — Cutaway view of the spherical sector model adopted for the circumstellar envelope. 

spherical sector model with R e — > oo has the peculiarity that the disk mass is infinite unless rj > 3.0. 
In what follows, the symbol rj will be retained for generality, but rj = 3.1 is used as a representative 
mean value in all the numerical computations. 



3. An Approximation for the Gray Polarization 

According to Fox (1991) except for a sign error in the first argument of the beta function 
B, the net polarization of light from a central star due to scattering of photons by free electrons 
(Thomson scattering) in a surrounding envelope may be roughly approximated for the specific case 
of the infinite spherical sector geometry as 



3(f? - 1) „ fv ~ 1 3\ 2 2 

Po = — B I — ^— , - 1 r e sin a cos asm i, 

where i is the angle of inclination of the rotation axis of the star to the line of sight and 
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is the total radial optical depth for electron scattering in the equatorial plane of the spherical sector 
in terms of the Thomson scattering cross section 
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Equation (2) is only a single-scattering approximation, based on the assumption that the disk 
is optically thin and neglecting multiple scattering. It includes a correction for the geometrical 
complication that the star is an extended source of light and cannot be treated as a point source 
(Cassinelli, Nordsieck, & Murison 1987, hereafter CNM), but it does not take into account the 
occultation of part of the scattering envelope by the star. 



So far the calculated polarization has no dependence on wavelength, since it comes from pure 
electron scattering. However, neutral hydrogen in the envelope absorbs some light both before 
and after scattering, and the envelope emission (which is treated here as unpolarized and is not 
considered to scatter) dilutes the polarization. It is therefore necessary to include these wavelength- 
dependent effects, since they combine to produce the slope of the continuum polarization and 
the well-known abrupt changes in the polarization at the H I ionization series limits, as first 
demonstrated by Capps, Coyne, & Dyck (1973) and later refined by McLean (1979). 

Assuming LTE ionization fractions and level populations, equation (5) gives the total wavelength- 
dependent opacity k(X) per gram of neutral hydrogen in the ground state (Aller 1963). Co is a fixed 
numerical factor, X(n) is the ionization energy from level n in units of kT e , and C se is the correction 
factor for stimulated emission. The three individual terms in the equation are a summation over 
the ionization edges of the first seven discrete energy levels, an integrated term for the combination 
of all the remaining levels, and a free- free absorption term (with Gaunt factors taken to be unity). 
See Appendix A for a discussion of NLTE corrections. 



4. Including the Wavelength Dependence 





where 
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If the free electron number density N e (r) and the electron temperature T e are known, then the 
number density N\(r) of neutral hydrogen atoms in the ground state can be found from the Saha 
equation: 

where 

The total radial optical depth for neutral hydrogen absorption in the equatorial plane of the spher- 
ical sector is therefore 



f' 00 ( B \ n 

r (A) = J K(\)m H N Q1 l-^j dr 

iVoimj/K(A)i?: 



^ (finite only for rj > 0.5) , (11) 

which acts as a wavelength-dependent attenuation factor to reduce the gray polarization. 
The volume emission coefficient of the envelope is 

f B \ 2r> 

j(\,r) = m H K(\)N 01 l-^j B(X), (12) 

where 

is the Planck function (not to be confused with the beta function in §3). The total luminosity of 
the envelope is then 



L(A) = 4vr / / / j(X,r)r 2 sin OdrdOdcf) 
Jo J%-a Jr* 

-N im H K(\)B(\)Bi sin a (finite only for 77 > 1.5) , (14) 



2r?-3 



- 7- 



which together with the stellar flux F*(A) produces a further wavelength dependence of the polar- 
ization. Theoretical values of the stellar flux F*(A) are taken from tabulated model atmospheres 
by Kurucz (1994). 

With the wavelength dependence included, our polarization estimate may now be written as 

PW = 1 + L(A)/(4^F,(A)) • (15) 

It should be emphasized that this equation is only an approximation (even if the disk is optically 
thin) , since it includes neither the contribution of scattering of the disk emission to the polarization 
nor attenuation of the direct starlight by the disk. Also the neutral hydrogen opacity is treated in 
a very crude way, using only the maximum radial optical depth in the equatorial plane. 

It may be of interest to calculate the total mass of the disk, which can be done approximately 
by simply counting pairs of electrons and protons under the assumption that it consists purely of 
fully ionized hydrogen: 



M = I / m H N 0e — r 2 sin 9drd8d<p 
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5. Synthesizing Broadband Data 

Adjustment of the envelope parameters in the theoretical model to give results consistent with 
observations might be expected to yield valuable information about the nature of Be disks. To 
compare with observations on the Johnson-Cousins UBVRI system, it is first necessary to convolve 
the theoretically estimated polarization p(A) from equation (15) with the bandpass characteristics 
of the optical system. 

As a first approximation, let us assume Gaussian filter transmission curves based on the stan- 
dard UBVRI wavelengths and fwhm (Bessell 1979), which are closely matched by the system 
actually used for the observations (see Table 1). 

Given the central wavelength A c , the standard deviation a can be calculated from the fwhm T 
since F = 2.354 a for a Gaussian distribution, so that the transmission function may be written as 



T(A) = exp 



A-A c 
2 Vr/2.354 



(17) 
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The theoretical total flux F(X) is due to the stellar component F*(X) plus the envelope com- 
ponent: 



F(\) = F m (\) + L(\)/(4irfi*). (18) 

Convolving the polarization in terms of flux with the transmission function, where \ m in an d 
\nax are, to a realistic approximation, the 10% response points, we find a theoretical expression 
for the expected polarization over the given bandpass: 

Ea-a x p(X)F(X)T(X)AX 
Etx min F(X)T(X)AX 



6. Adjustable Parameters & Model Fits 

With the model in hand, an interactive graphical computer interface was designed to allow ad- 
justing the input parameters by trial and error to fit UBVRI polarization measurements (McDavid 
1999) of eight Be stars. Solutions obtained this way should be viewed with some caution since 
they may not be unique: different geometrical distributions of scatterers can produce the same net 
polarization. Moreover, given the gross approximations involved, only order-of-magnitude results 
should be expected. 

The model disk has three adjustable parameters: the maximum electron number density iVoe, 
the angle of inclination i of the rotation axis to the line of sight, and the opening half-angle 
a. Additionally required fixed parameters for each individual star are the radius i?*, the effective 
temperature T*, and the surface gravity log g, which were estimated from Table 1 of Collins, Truax, 
& Cranmer (1991) based on spectral types from Slettebak (1982). The electron temperature of the 
disk was then fixed at T e = 0.75 T* and the appropriate Kurucz flux table chosen to match T* and 



Table 1. Filter System Parameters 



Filter 


Effective Wavelength 


Bandpass (fwhm) 




(A) 


(A) 


U 


3650 


700 


B 


4400 


1000 


V 


5500 


900 


R 


6400 


1500 


I 


7900 


1500 
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log g- 

The fixed disk temperature T e affects the wavelength dependence of the polarization through 
its influence on the neutral hydrogen opacity by setting the degree of ionization and the populations 
of the excited states (eq. [5]). It influences not only the hydrogen absorption optical depth (eq. [11]), 
but also the disk luminosity (eq. [14]), which is another contributor to the wavelength dependence 
of the polarization. As a result, T e is important in determining the polarization Balmer jump 
and the slope of the polarization over the Paschen continuum, both of which generally increase at 
lower temperatures. Appendix A explains the NLTE corrections to the level populations which are 
necessary to fit these features for some of the program stars. 

It is instructive to summarize individually the effects of the three adjustable parameters on 
the behavior of p{\) according to equation (15). 

(1) The overall degree of gray polarization is directly proportional to N £ q through the factor 
r e in equation (2). However, simply increasing N e Q does not always result in an overall polarization 
increase, because it also increases the attenuation and adds to the wavelength dependence of the 
polarization through r a (A) (eq. [11]) and L(X) (eq. [14]). 

(2) The earliest basic electron scattering models (e.g. Brown & McLean 1977) gave a sin 2 i 
dependence of the polarization on i, and further refinements have not qualitatively changed the 
simple picture of Figure 1 in which a face-on disk shows no net polarization because of symmetrical 
cancellation, while the maximum asymmetry of an edge-on disk results in the maximum polariza- 
tion. It may seem likely that i is poorly determined in this model since the same polarization might 
be obtained, for example, by decreasing i while increasing a, keeping iV e o constant. Experiments 
do not bear this out, though, because a thicker disk has a higher luminosity, which changes the 
wavelength dependence of the polarization enough to ruin the fit (eq. [15]). Another helpful con- 
straint on i is the presence or absence of shell lines, which are thought to be absorption features 
in the spectral line profiles due to circumstellar material in the line of sight for nearly edge-on 
equatorial disks. Since Be stars are generally rapid rotators, spectroscopic measurement of the 
projected rotational velocity v sin i can also suggest roughly the value of i. 

(3) Figure 3 shows how the model polarization p B in the B passband varies as a increases from 
0° to 90° using a typical set of parameters. The polarization first rises as the widening disk presents 
an increasing number of scatterers, then begins to decline due to the buildup of enough density 
at high latitudes to cancel the polarization from the equatorial regions. The model of Waters 
&: Marlborough (1992) also shows this behavior, which becomes of practical concern because the 
hump-shaped graph allows for both a "thick disk" and a "thin disk" solution, producing identical 
polarization at a values on opposite sides of the peak. For equal densities, however, the thick disk 
will have a higher luminosity than the thin disk, giving the polarization a recognizably different 
wavelength dependence. 

Figures 4-11 present possible (but not necessarily unique) model fits for the eight program Be 
stars, with the parameters summarized in Table 2. The observational data points, plotted as open 
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Fig. 3. — The dependence of the model blue polarization p B on the opening half-angle a of the disk using a typical 
set of parameters. 

squares with vertical error bars, are intrinsic polarization, corrected for the interstellar component 
by Stokes vector subtraction (McDavid 1999). The corresponding model result is plotted as a solid 
line, using open circles for the broadband values with horizontal line segments to show the filter 
passbands. Several additional quantities derived from the model are also included in Table 2 for 
each fit: the maximum polarization pmax and its wavelength X(pmax), the maximum of the H I 
absorption optical depth T a max and its wavelength \{T a max), the electron scattering optical depth 
r e , the total envelope mass M, and the NLTE departure coefficients 62 and 63 for the populations 
of the first two excited states of the neutral hydrogen in the disk. Refer to Appendix A for a 
discussion of NLTE considerations. 

7. Test of Accuracy 

Since the model presented here is such a crude approximation, it is important to compare 
its results with those from a more sophisticated model to investigate, at least qualitatively, the 
accuracy we can expect. The Monte Carlo code of WBB is one such standard, and it is based on 
the same set of input parameters except that the density input is the electron scattering optical 
depth rather than the electron number density, which requires only a minor conversion. Ideally 
we would like to determine WBB fits to the program star observations, and then quantitatively 
evaluate the simplified model by direct comparison of the fit parameters. Unfortunately, because of 
the slow convergence of the Monte Carlo method, deriving a WBB fit is a complex process requiring 
substantial amounts of cpu time even on a fast computer. Otherwise there would be no reason for 
interest in the kind of quick approximation that is the subject of this paper. 
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Fig. 4. — Fit of a spherical sector disk to the UBVRI polarization of 7 Cas. 
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Fig. 5. — Fit of a spherical sector disk to the UBVRI polarization of (f> Per. 
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Fig. 6. — Fit of a spherical sector disk to the UBVRI polarization of 48 Per. 




2000 4000 6000 
A (A) 



8000 10000 



Fig. 7. — Fit of a spherical sector disk to the UBVRI polarization of £ Tau. 
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Fig. 9. — Fit of a spherical sector disk to the UBVRI polarization of \ Oph. 



-14- 



2.0 
1.6 



1.2 

0.8 

0.4 
0.0 

2000 4000 6000 8000 10000 

A (A) 

Fig. 10. — Fit of a spherical sector disk to the UBVRI polarization of n Aqr. 
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Fig. 11. — Fit of a spherical sector disk to the UBVRI polarization of o And. 
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Table 2. Be Star Modeling Data 





7 Cas 


4> Per 


48 Per 


C Tau 


48 Lib 


X Oph 


7r Aqr 


o And 


Spectral Type 


B0.5 IVe 


B1.5 (V:)e-shell 


B4 Ve 


Bl IVe-shell 


B3:IV:e-shell 


B1.5 Ve 


Bl III-IVc 


B6 III 


Date of Obs. 


1997.96 


1997.96 


1997.96 


1994.94 


1998.55 


1995.49 


1991.59 


1998.55 


R*(Rq) 


6.5 


6.0 


4.5 


7.0 


6.0 


6.0 


9.0 


5.5 


n (K) 


25,000 


25,000 


17,000 


25,000 


20,000 


25,000 


25,000 


14,000 


log g 


4.0 


4.0 


4.0 


4.0 


3.5 


4.0 


3.5 


3.5 


-*-r- / J rt 1 O Q \ 

7V 0e (10 12 cm- 3 ) 


9.00 


5.00 


3.00 


4.20 


4.60 


6.20 


6.50 


2.00 


T e (K) 


18,750 


18,750 


12,750 


18,750 


15,000 


18,750 


18,750 


10,500 




2.0 


10.0 


3.0 


6.0 


4.0 


3.0 


3.5 


5.5 


i(deg) 


52.0 


80.0 


55.0 


80.0 


80.0 


41.0 


80.0 


80.0 


pmax(%) 


0.60 


2.41 


0.25 


1.51 


0.94 


0.42 


1.66 


0.54 


X(pmax) (A) 


3661 


3862 


3862 


3661 


3812 


3661 


3661 


3862 




0.94 


0.27 


0.23 


0.22 


0.44 


0.41 


0.68 


0.22 


X(T a max) (A) 


3636 


3636 


3636 


3636 


3636 


3636 


3636 


3636 


T e 


1.29 


0.66 


0.30 


0.65 


0.61 


0.82 


1.29 


0.24 


M(1O- 9 M ) 


3.07 


6.69 


0.51 


5.37 


2.47 


2.50 


10.3 


1.14 




1.69 


1.69 


1.05 


1.69 


1.35 


1.69 


1.69 


0.70 


&3 


0.29 


0.29 


0.23 


0.29 


0.26 


0.29 


0.29 


0.19 
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For the time being, though, we can at least make use of the existing WBB fit to £ Tau. Figure 12 
shows examples of p(X) for the quick approximation (solid line with broadband filter points marked 
by open circles) compared to the full optically thick multiple scattering WBB simulation (dashed 
line connecting filled circles). Errors in the WBB data points are on the order of 0.03%, limited 
by computing time. The observed UBVRI data are plotted as open squares with error bars. Using 
the parameters R* = 5.5 R & , % = 19, 000 K, T e = 15, 000 K, a = 2?5, and i = 82?0 from the WBB 
solution as input to both models, the six plots show what happens as r e is increased by steps from 
0.10 to 3.00 (the best fit value according to WBB). 

This test shows that the optically thin approximation agrees quite well with the WBB calcula- 
tion up to and including an optical depth of 1.00. It is also encouraging that for £ Tau the best fit 
of the approximate method as shown in Figure 7 and Table 2 has iV e o = 4.20 x 10 12 cm" 3 , which 
is well within the range of 2.38 x 10 12 < N 0e < 7.54 x 10 12 cm" 3 found by Waters, Cote, & Lamers 
(1987) in their infrared study. 

However, it must be noted that there is a possibly serious disagreement between the WBB 
model and the simplified model presented here with regard to the best-fit values of r e and a for 
C Tau. While the WBB disk (r e = 3.00 and a = 2?5) is optically thick but geometrically thin, the 
simplified model disk (r e = 0.65 and a = 6?0) is optically thin and geometrically thicker (although 
it still qualifies as "geometrically thin" ) . 

In Figure 12 the approximate model fails to produce a polarization large enough to fit the ob- 
servations of C Tau by increasing r e at constant a because it becomes dominated by H I opacity. As 
shown in Table 2, this problem is seemingly avoided by increasing a instead of r e , but the resulting 
optically thin solution could possibly be ruled out on other grounds, such as the associated infrared 
emission as discussed by WBB. This is only one cautionary example emphasizing the importance 
of independent observational consistency checks in the interpretation of Be-star polarimetry data. 

8. Conclusion 

In summary, the model presented here appears to give reasonable results for optically thin 
cases, but we should only trust it for order-of-magnitude accuracy because it is based on so many 
simplifying assumptions. Even so, it may be useful for estimating trial values of the parameters as 
starting points for rigorous fitting with more elaborate models. 

An interesting result from applying the model to observations of the eight program stars is 
that all of them may be fit with geometrically thin disks, with opening half-angles of ten degrees or 
less. This adds tentative support to the statistics in favor of thin disks as presented by Bjorkman 
& Cassinelli (1990) and later discussed at length by WBB. If most Be disks prove to be similarly 
flat, it will surely have important implications concerning their formation process. 

It may seem discouraging that this model will never find optically thick solutions even if they 
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Fig. 12. — A comparison of the simple approximation (solid line with open circles) against the WBB solution 
(dashed line with filled circles) for a range of electron scattering optical depths r e . The observational data points for 
C Tau are shown as open squares. 
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physically exist. And with three adjustable parameters, there might be any number of apparently 
valid optically thin fits to a set of UBVRI data points, rendering the entire exercise inconclusive. 
Nevertheless, there are good reasons for a more optimistic point of view. Due to the efforts of 
many investigators, as documented, for example, in the references in §1 of this paper, the range of 
physically allowable parameters has been narrowed significantly. And within even the most liberal 
of these constraints, experience with a quick, interactive trial-and-error system based on even the 
elementary physics presented here will show that it is surprisingly difficult to find multiple solutions 
with qualitatively different sets of parameters. I will be glad to share my IDL modeling program 
with interested users on request. 

Thanks to Joe Cassinelli for critical reading and commentary on early drafts of this paper, 
and especially for suggesting a method to estimate the NLTE departure coefficients for the neutral 
hydrogen in a Be disk. Kenny Wood and Barbara Whitney also provided very helpful discussions, 
and Jon Bjorkman kindly gave me a version of the Monte Carlo polarization code. I especially 
thank the anonymous referee for pointing out many serious errors in the original manuscript and 
for patiently explaining the necessary corrections. 



A. Estimating the NLTE Departure Coefficients 

Fits to the polarization Balmer jump were difficult and in some cases impossible without NLTE 
corrections to the level populations of the first two excited states of H I. According to the com- 
bined Saha and Boltzmann equations, the maximum population of the first excited state normally 
occurs at a temperature of about 10,000 K and declines rapidly toward higher temperatures with 
increasing excitation and ionization. At temperatures closer to 20,000 K appropriate for Be disks, 
the LTE value of the absorption coefficient at the [/-band is sometimes too small to match the large 
polarization Balmer jumps which are commonly observed. NLTE correction is applied only to the 
first two excited states because they dominate the opacity at optical wavelengths for the conditions 
of density and temperature of interest. 

Since only averages over a line of sight are required for this very approximate polarization 
model, departure coefficients were simply estimated from the calculations of CNM. Their Figure 3 
shows mean population parameters q n for n = 2 and n = 3 defined by 



QnNLTE 



f£q n (r,T.)N?dr 

(A1) 



based on the detailed NLTE analysis of a stellar wind immersed in the radiation field of a central 
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star, with statistical balance between photoionization and radiative recombination. The departure 
coefficient for level n is then the ratio of the NLTE and LTE values of q n . 

The calculation of q nLTE is done as follows: 



/ N n {r)dr = n 2 N 01 e x W- x W (^) dr 

Jr* Jr, \ r J 



n N 01 R* x(n)-x(i) (A2) 
27] -I v ' 



and 



f-co f-oo / tj \ 2r; 



NlR* 
27] -I ' 



(A3) 



so that 



QnLTE 



= ^0i n 2 e X(n)-X(l) 

= (2TTm e kT e )- 1 - 5 h 3 n 2 e x ^ . (A4) 

The population parameters Q2NLT E an d Q3NLTE ma Y De estimated by extrapolation of power 
law fits to the disk temperatures (10,000-20,000 K) and mass-loss rates (~ 10 -8 M Q yr~ l ) char- 
acteristic of Be stars, based on Figure 3 of CNM: 

g mTE = i.iixio- 8 r; 2 - 83 , (A5) 



Wt£ = 4.59x10- 13 T- 2 - 02 . (A6) 

The approximate departure coefficients are then 62 = Q2NLTE/Q2LTE an d ^3 = Qsnlte/Qslte-, 
so that the corrected level populations are roughly N2N lt e = b2 x N2LTE and N^nlte = b3 x 
N3lte- These correction factors are applied to the appropriate terms in the summation inside the 
parentheses in equation (5). The corrected form of «(A) is used only in calculating the absorption 
optical depth and not in calculating the emission coefficient, since recombination and free- free 
emission are LTE processes. 
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